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CJ '. Abstract 

I ■ We perform a minisuperspace analysis of an information-theoretic 

bi) nonhnear Wheeler-deWitt (WDW) equation for de Sitter universes. 

The nonhnear WDW equation, which is in the form of a difference- 
CN ■ differential equation, is transformed into a pure difference equation for 

the probability density by using the current conservation constraint. 



> 

T-|- . In the present study we observe some new features not seen in our 

C^ ■ previous approximate investigation, such as a nonzero minimum and 

maximum allowable size to the quantum universe: An examination of 
t^ , the effective classical dynamics supports the interpretation of a bounc- 

ing universe. The studied model suggests implications for the early 
universe, and plausibly also for the future of an ongoing accelerating 
phase of the universe. 

X: 



1 Introduction 

While quantum physics is expected to play an important role near the Big 
Bang when the universe was small, there are a number of differing opinions 
on how to quantise the classical universe, see reviews in [H [2], leading to 
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different approaches to quantum cosmology. One issue of interest in quantum 
cosmology is whether classical singularities are resolved [3]. In our previous 
study [5! we found the answer to be affirmative in the information-theoretic 
approach, as is the case also in many other approaches [U 121 13]. 

The philosophy of the information-theoretic approach, also known as the 
"maximum uncertainty (entropy) method" , is that one should minimise any 
bias when choosing probability distributions, while still satisfying relevant 
constraints [5j. While this approach originated in statistical mechanics, it has 
wider applicability such as motivating the structure of the usual Schrodinger 
equation [6] and its potential nonlinear modifications to model new short 
distance effects[71 EJ E]. 

In our previous study [1|, which we briefly review in the next section, we 
treated the nonlinearity of the modified WDW equation ([7]) perturbatively 
and studied a truncated linearised equation with an effective potential. In 
this paper we present a full investigation for the case of a de Sitter universe. 
The de Sitter universe is a reasonable model for early stages of inflation and 
furthermore, if the current accelerating phase of our universe continues, then 
at late times it can again be well approximated as a de Sitter universe. 

Treating the universe as an isolated system, one may apply quantum 
mechanics to the whole de Sitter universe regardless of its size; in other 
words, study the wavefunction of the universe. Our main intent in this paper 
is to see how the classical de Sitter dynamics is modified by information- 
theoretically motivated corrections to the WDW equation, extending the 
previous [4J approximate investigation. 

Consider the Einstein- Hilbert action for a FRW universe with the cosmo- 
logical constant A = S/oq modeling inflationary sources in the early universe. 



S= I dtL = - I dtN 
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where < a(t) < oo is the dimensionlesci scale factor, N the lapse function, 
fc = 0, ±1, and we have taken h = c = 1. Varying with respect to A^ and 
choosing the A^ = 1 gauge gives us the Friedmann equation 

d^ + (/c - a^) = . (2) 

We have set oq = 1 as the results for other values can be obtained by scaling. 



^As in [4;, the physical scale factor aphys — Pa, where /3 oc Ip, the Planck length. 



The flat, k = 0, geometry has the expanding classical solution 

a = exp(t) (3) 

which implies an arbitrarily small universe, a — ?• 0, at early times t — )• — oo. 
Expected quantum effects are studied in standard minisuperspace quantisa- 
tion by promoting the canonical momentum 11 = ^ = —da to an operator, 

n -)■ n = -i^, leading to the WDW equation [10], 



'da ' 



d^ 



da^ 
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(4) 



The solution of (j4]) representing an expanding universe is given by the 
Hankel function 
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as a — )■ oo , 
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as one can verify by noting that the negative value for the momentum 11 = 
—ad of the asymptotic form corresponds to d > 0. Notice that the , solution 
([S]) of the linear WDW gives a nonzero probability for a universe of size a = 0, 
and hence for a Big Bang. In a previous perturbative study [1] we showed 
that a modified WDW dynamics screened the a = region, leading to the 
quantum creation of a universe through tunneling [10] . 

We postulated in Ref.[3] that if there is some new physics at quantum 
gravity scales, it may be modeled, within the information theory framework, 
by a modified WDW equation |7] which can be derived by extremising the 
usual WDW Lagrangian while also maximising the relevant uncertainty mea- 
sure. We state and explain first the equation before motivating it further 
below: 

-a^ + F{p) 



da^ 
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F{p) = Qnl-Q, 


with 
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and 

^ ^, da? ' ' ' 

Here p(a) = ip*{a)%lj{a) is the probability density and p±(a) = n{a±C), where 
the dimensionless parameter ^ > is the nonhnearity scalqj. The hnear 
theory is recovered as C — ^ while the other parameter < r/ < 1 labels a 
family of nonlinearisations. Notice that the modified WDW equation is still 
invariant under a scaling of the wavefunction, ip — )■ \ip, and so the solutions 
of the equation do not depend on the normalisation of the probability density. 
The appearance of p± in ([7]) implies a nonlocality of the modified WDW 
equation which takes the form of a difference-differential equation. While 
the equation itself looks complicated, it follows from a relatively simple La- 
grangian involving the Kullback-Leibler (KL) information measure [71 H]: 
The nonlinear piece F in Eq. ([7]) arises from maximising, in the spirit of the 
maximum uncertainty method, the following term in the action 

/pid) 
P(o)ln— — —r^da. (11) 

The Kullback-Leibler measure is clearly a relative uncertainty measure 
which generalises the usual Gibbs-Shannon entropy of statistical mechanics, 
and which reduces, as C ~^" 0, to the "Fisher" information measure used in 
Ref . [6] for motivating the usual Schrodinger equation within the information- 
theoretic framework. Indeed, as discussed in [7], the (77- regularised) KL mea- 
sure is probably the simplest nonsingular measure which interpolates between 
the Shannon and Fisher measures and which keeps some desirable properties 
such as the scale invariance mentioned above. 

In physical terms, the nonhnearity scale C, may be interpreted as a kine- 
matic implementation of the resolution at which the coordinates become 
distinguishable \l\. Clearly for C 7^ the dynamics given by Eq.Q will be 
modified from the usual case. 

The nonhnearity in ([7]) was originally explored in quantum mechanical 
systems [7J, and its various perturbative and non-perturbative properties 
studied in Refs.[7l El [9l [H]. However as experimental constraints simply 
place limits on the size of the nonhnearity for simple quantum mechanical 
systems [7], we then applied the modified quantum equation to cosmology in 
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Ref . [1] and found encouraging results for singularity resolution. In this paper 
we continue our study of the consequences of the information-theoretically 
motivated nonlinear WDW equation (jTj). 

The outline of this paper is as follows. In the next section we review 
the perturbative results of Ref. |1] on how the nonlinearly corrected quantum 
dynamics can avoid the a = possibility that is present in Eq.([5]). Then in 
Section(3) we show how to transform ([7]) into a purely difference equation 
for p which is solved numerically in Section(4) and studied analytically in 
Section(5). In Section(6) we discuss the effective classical dynamics suggested 
by the nonlinear WDW equation and elucidate the physical meaning of the 
nodes of the wavefunction. In Section(7) we discuss the A = case and 
conclude in Section(8). 

2 Review of perturbative treatment 

If the nonlinearity F is weak, it may be expanded perturbatively for C ^ 1, 
giving to lowest order 

F(p) = C(3-4r^)/(p) + 0(a, (12) 

where 

m = ^,{2p" - 3p"p) . (13) 

In this approximation Eq. ([7]) is 



da^ 



a4 + C(3-4r7)/(a) 



i){a) = (14) 



which we may solve by iterating about the unperturbed solution ([5]): At low- 
est order one calculates po = ipQ^po and then /o(a), giving a linear Schrodinger 
equation with an effective potential 

K// = -a' + C(3-4r/)/o(a). (15) 

The perturbative approximation (1121) around the linear solution requires 
not just C ^ 1 but also that /o(a) be slowly varying, which is indeed the case 
since for a — )■ oo we have /o ~ 1/a^. In particular, there are no singularities 
from nodes of the wavefunction in the expansion of F{p) around po, unlike 
the case of quantum mechanical systems studied in [S] . 



As shown earlier in [1], for 77 < 3/4 the nonhnearity fo{a) forms an 
effective potential barrier, a finite size universe coming into being through 
quantum tunneling \IU\. In other words, in the modified classical dynamics 
a backward evolving classical universe will experience a bounce instead of 
shrinking to zero size. 

In the sections below we venture beyond the perturbative analysis of 
Ref . [4j to uncover some new features of Eq. ([7]) . 

3 Exact difference equation for p 

Writing the exact wavefunction in terms of its amplitude and phase, ip = 
y^e**^, with p and 5" real, the imaginary part of the WDW equation ((71) is 
then the continuity equation 

which can be solved to give 

"fa— (1^) 

The constant current a is fixed by requiring our nonperturbative solution 
approach the asymptotic forno of the solution for the linear theory (j6]) near 
some large a, which corresponds, for example, to "now" . This gives 

a = ^. (18) 

71 

Next, we use ( ITT]) in the real part of the nonlinear WDW equation ([7]) to 
eliminate the derivatives of S, giving a purely difference equation for the 
probability density: 



a' - Qnl (19) 
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""However as we shall see, for even larger a, beyond the matching point, the nonlinear 
WDW solution will eventually deviate from the form (j6]) indicating an eventual departure 
from the current classical evolution. 



In the derivation we factored a common y^ from both sides of (I19p : 
as discussed at the end of Sect. (5.1), this does not affect the results even as 
p — )■ 0. The difference equation f lT9|) relates the adjacent values of probability 
density p_,p and p+ which are separated by the step size (, the nonlinearity 
scale 

It is important to note that the variable a is still continuous, as is p{a). 
It is just that equation ( 1T9|) places non-local constraints on the p{a) values. 
Thus although the solutions to be discussed below are on a lattice of step 
size C as determined by ( 1T9|) . it is to be understood that the region between 
the discrete set of points is continuously connected. 

The equation ( IT9l) for the probability density does not depend on the sign 
of a and hence describes both possibilities, either an expanding or contracting 
universe; the specific wavefunction describing either possibility does depend 
on a. 

As C — ^ 0, —Qnl becomes —Q > 0, Eq.f jTOj) . which is positive definite for 
the linear WDW equation solution ([5]) and hence the right-hand side of ( !T9|l 
becomes positive definite, just as the left-hand side already is. However, for 
C 7^ the difference equation ( IT9|) restricts the range of a as we are required, 
by definition, to preserve the positivity of the probability density. When 
starting with p{a) > at some initial point and using the difference equation 
to move forward or backward, it is possible that one reaches a point at which 
p = 0, beyond which p < or becomes complex, as the equation ( 1T9l) by itself 
does not guarantee a real or positive p±: the kinematic constraint p > must 
be self-consistently imposed. 

In Sects. (4,5) we interpret the occurrence of p = as delimiting the range 
of allowed a values and hence on the size of the universe: The consistency 
of such an interpretation will be seen by examining the effective classical 
dynamics in Sect. (6). 

In summary, the kinematic physical constraint p > 0, when imposed on 
the dynamical difference equation (TT9|) . constrains the size of the universe. 
(Such a constraint does not occur for usual quantum mechanical systems, see 
later) . 

4 Numerical Results 

The difference equation (TT9|) is easily solved by specifying two initial values 
which we fix by requiring the large universe to be close to the asymptotic 



form (jH]) given by the usual linear WDW equation near a = 5. Since Eq. fll9p 
easily gives p_ explicitly in terms of p and p+ , hence if the later two variables 
are initially fixed then a direct backward evolution of the equation gives the 
values of p for smaller a. The numerical accuracy was set at 16 figures which 
allows the features of p discussed below to be unambiguously distinguished. 

However in Eq. flT^ p^ cannot be written explicitly in terms of the other 
two values, so values of the probability density forward from the starting 
point were obtained by solving the implicit equation using Newton's method. 
We checked the accuracy of Newton's method by making several consistency 
comparisons; for example, using the end points of the forward steps as start- 
ing points for the direct backward evolution and comparing the two curves. 
Identical results were obtained using the "Solve" function in MATLAB. 

We summarise below the key results. Unless otherwise stated, rj = 0.5. 
For the nonlinearity scale, we explored the range 10~^ < ( < 1 though not 
with the same degree of detail for every feature. 

1. For low nonlinearity, C — > 0+, starting from a = 5 and moving back- 
wards towards a = 0, p remains positive and close to the p for the 
linear theory. These results, including the effective potential Vejf = 
—a^ + F{p) which develops a barrier, Fig.(l), agree with those obtained 
using the lowest order perturbation theory in Ref. [4]; the barrier being 
smaller for smaller (. (For the evaluation of the effective potential us- 
ing the discrete p data we used the central difference approximation for 
the second derivative in Q f lTU]) . which is sufficient for the low ( values 
encountered.) 

However, moving forward from the starting point, the probability den- 
sity will eventually become zero at somqj finite value which we label as 
dmax- The consistency of interpreting amax (and amin below) as their 
labels suggest will be seen in Sect. (6). The value a^ax depends on r], (, 
for example amax = 14.4 for rj = 1/2 and ( = 0.005 for the initial 
conditions used; Figures (2,3) show the curves. The general trend is 
that amax decreases with increasing ( ioi ( < (c (as defined below). 

Thus, unlike the perturbative study in Ref.[4], the non-perturbative 
results show a maximum allowable size to the quantum universe. For 
low values of (, amax fits a power law in (, Fig. (4). 



^In the numerical work we locate the first point where p becomes negative or complex. 
Due to the lattice nature of ([T9|) . the actual point where the probability density first 
vanishes will be between two lattice points. 



2. As the nonlinearity ( increases beyond a critical value (c, P becomes 
zero at some aminijlX) during the backward evolution, in addition to 
vanishing at a forward point amaxivX) during the forward evolution. 
For example, Cc = 0.049 for rj = 0.5 and the initial conditions used. 
Fig. (5). In the range C ^ Cc the trend is that increasing ( leads to 
increasing a^ax and decreasing amin- 

The existence of amin meanqj that the quantum universe has a nonzero 
minimal allowable size: This is the second new result not seen in the 
perturbative treatment of Ref . |H] . Such a possibility, of a minimum size 
to space determined by the quantum nonlinearity F, was earlier noted 
inRef.[9j. 

3. As ( increases from zero, oscillations become apparent in the probabil- 
ity density at larger values of a, in contrast to the probability density 
for the linear theory ([5]) which is monotonically decreasing. These os- 
cillations are discussed in more detail in the next Section. 

4. We checked that the new features observed in our numerical study of 
flT9|) are robust to changes in initial values for p or the starting point 
a used. For example, we solved the equation for starting points near 
a = 10 and used starting values for two of the adjacent p's in fITIJl) to be 
the same small value; the results again show the existence of amin, amax 
and the oscillations though the specific numerical values differ. 

5. All the above results are qualitatively the same for a lower value rj = 0.2 
that we checked. However for larger t] = 0.9, though there is still an 
amax and a Cc there is no potential barrier at small ( in agreement 
with Ref.[l] which showed the barrier to be present at low ( only for 
T] < 0.75. 



5 Analysis 

The features observed in the numerical study may be understood through 
various analytical approximations which we discuss in this section. 



^In some cases, a continued evolution beyond amin or Umax, where p is negative or 
complex, leads to new regions where p becomes positive again. We assume that the 
wavefunction must be continuous and so display only the positive component of p around 
the initial point. 
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5.1 Existence of 0^,17 and a 



max 



As an illustration, a relation for amax can be estimated for ( small as follows: 
Set p-|_ = at a = amax in f lT9|) and as C — ^ estimate ip = ip{amax — C) ~ 
—(i/j' = — 7C so that p ~ (tC)^- Similarly p_ ^ (27C)^ where 7 is the slopqfl 
of the wavefunction at a. Then 

Note that the slope 7 is in general a function of ( and ctmax but if that 
dependence does not overcome the explicit ( factors in fl20p then we see that 
o-max — )■ cxD as C — ^ 0. The numerical results for amax are shown in Fig. (4). 
For ( small we find amax oc -j=. 

While setting p+ or p_ to zero in Eq. (IT9l) is unproblematic, it appears that 
p — )■ leads to a divergence through the ln(p) term. However re-arranging 
f fT9|) shows that p — )■ leads to the following consistency relation 



2ri'^C^a^ = limp2 



?7P+ r^p. 



l-77)p + 77P+ {l-r])p_ + rip 



(21) 



which implies that the expression in square-brackets must develop a — 1/p^ 
divergence through either p^ or p_ becoming negative. That is, as p — )■ 0, 
one of the adjacent points enters the unphysical region beyond either amax 
or amin- As discussed earlier, this may be interpreted as implying that the 
quantum universe within this nonlinear WDW framework is boundedjj. 

5.2 An Exact Solution 

An exact analytic solution of flTOl) . which illustrates the occurrence of both 
an amin > and amax, is obtained by taking three lattice points with p_ = 
p+ = and p = p{am) > at the mid-point am- Then flT^ gives p in terms 

oiamX^V- 

^^ 4 1 1 

p2 rn 2ri^Q 1-1] ^ ^ 



®We are assuming here that it is the wavefunction which is smooth near the node. If 
instead the probabihty density is smooth then an expression different from but similar to 
Ea. (|20)) can be derived. 

^The situation is different for the analogous difference equation for bounded potentials 
in quantum systems [II HI HI [11] : In those cases such a constraint does not arise. 
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Since the right-hand side of (1221) must be non-negative, this sets one con- 
straint. Also since a^ — C = cimin > this gives the second constraint 
dm = C + O'min- Thus fl22|) imphcs 



C' (C + dminY > A In -^ . (23) 

17]'^ 1 — 7] 

One can have large universes, up to size a^m + 2C, by taking 77 — > 0+ or 1~. 

5.3 Oscillations in the Probability Density 

Return to the nonlinear WDW equation ([7]) and now write the the wavefunc- 
tion as exp {6/2 + iS) for real 6, S and treat the nonlinearity F as one piece 
rather than separating Qnl from Q. As in Sect. (3), one may eliminate S 
and obtain the following equation, 

29" + {e'f + 4(a^ - a^e-^') + e = , (24) 

where e represents the terms from the nonlinearity F which we treat in the 
analysis below as small perturbations to the original linear WDW equation; 
the prime means d/da. Then, if 9q is the solution to Eq. fl24|) for e = 0, the 
full solution may be written 6 = 6q + 5 and this substitution in ( 121|) gives, to 
lowest order in e and 5, the equation for the fluctuations 

6" + e'^6' + AaHe-^'^' = . (25) 

Since the probability density po corresponding to the solution of the linear 
WDW (jlj) is monotonically decreasing, so 9^ < Q and hence the fluctuation 
equation ( l25l) describes oscillations which increase with a (anti-damping). As 
Po '^ 1/a? for large a, therefore 9q oc — 1/a and the anti-damping eventually 
vanishes. Furthermore, the last term of f l25|) implies that the wavelength of 
oscillations decrease with a. 

This analysis explains why in our numerical solutions the oscillations were 
not visible at low a: their initial amplitude was small and the wavelength 
large. The oscillations became manifest only as their amplitude increased and 
wavelength decreased at larger although the amplitude eventually stabilises, 
the wavelength keeps decreasing^, see Fig. (6). 



^^Notice that the leading order analysis does not refer to the specific form of F. Any 
small perturbation of the original linear WDW equation is expected to give rise to quali- 
tatively similar oscillations in the probability density. 
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Including sub-leading terms in (12^ would give the oscillations a depen- 
dence on the parameters 77, C as we have also observed numerically. As 5 — )■ 
when (^ — )■ 0, the amplitude of the oscillations must also be of order ( or 
smaller. The ratio p/po = exp(5) ~ 1 + 5 then implies small oscillations 
about the value 1 as illustrated in the actual numerical results of Fig. (5). 

We have also studied the truncated differential equation flT^ obtained 
by expanding F to lowest order in (. Direct solution of that nonlinear dif- 
ferential equation shows oscillations for p, as does the first order iterative 
version described by the linear Schrodinger equation with effective potential 
(fTSjl . However neither of these approximations to Eq. (fT9|) display the limiting 
values amax,cimin sccu in the nonperturbative difference equation (TT^ . 

The oscillations we have observed here do not seem to be related to the 
oscillations of the non-perturbative mode of a similar nonlinear Schrodinger 
equation studied in Ref. |l7J: the latter oscillations were of constant wave- 
length. 

6 Effective Classical Dynamics 

As the WDW equation is time-independent, it cannot describe the time- 
evolution of a universe. Furthermore the probability density refers to the 
likelihood of observing a universe of a particular size in an ensemble. There- 
fore, in order to gain some insight into the dynamics of a single universe, in 
particular near amax and flmm, we return to the classical domain. Near the 
nodes, we see from fITTj) that dS/da is large and as argued in Refs. [121 ^Q\ 
it may be identified with the classical momentum, as in the discussion after 
Eq.®. 

In that approximation the effective classical dynamics corresponding to 
([7]) is described by the modified Friedmann equation 

a^d^ + Veff = (26) 

where V^jf = —a^ + F{j)) is the effective potential and p is the solution of 
the nonlinear WDW equation. Note that the modified Friedmann equation 
is again a difference-differential equation. In Ref.jl] we treated fl26|) pertur- 
batively by expanding F to lowest order in C, but in the discussion below we 
keep the full form of F. 

If the nonlinearity F is weak and non-singular then it cannot overcome 
the classical piece — a^, especially at large a. However near a = a^, > 0, where 
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the wavefunction vanishes, p ~ (a — a*)^ and so from Sect. (5.1) we see that 
Qnl ~ ~l/p^ near a*. Similarly, from Eq. inU]) . Q is also likely to be large 
and diverging, thus making F{p) = Qnl — Q large and possibly positive. 
(The enhancement of the nonlinearity near nodes of p was earlier noted in 
Ref.|8]). 

A numerical examination of V^ff shows that between amin and amax it 
is real and negative. In some cases Vg// remains real but becomes positive 
as one approaches either amin or amax thus forming a potential barrier there 
which imphes, through ( 126|) . a bounce: For example, approaching amax from 
the left, and assuming a Taylor expansion of Eq. (!26|) near that point, one 
obtains (dropping positive constants) 



V(^max — « , (27) 



"max- 



showing that a = and a < as a — )■ a^ 

However for quantum solutions which have both an amin and amax, we 
found that if there is a potential barrier in the effective classical dynamics 
at one end, say near amin, then the effective potential at the other end, near 
amax is complex. The complex Veff implies that p has reached an unphysical 
value (negative or complex). Since Qnl is defined only at a discrete set of 
points, it is possible that F and hence Veff is not analytic near a = a^,. For 
example, the \np ~ ln(a — a*)^ term in Qnl might be dominant leading to a 
transition from a negative to a complex, and hence unphysical V^ff. In such 
situations the approximation ( 126|) has broken down near that point. 

In Ref . [H] , in the context of a different quantum cosmological model, the 
authors found that locations of singularities in the classical dynamics were 
where the wavefunction vanished in the quantum dynamics. In some sense 
we have found a complementary situation, whereby nodes in the probabil- 
ity density of a nonlinear quantisation implies, in some cases, an effective 
classical dynamics with singularities near the nodes. 

However it is possible to observe a cyclic universe in the effective classical 
dynamics if we take C < Cc and rj < 3/4, for which there is no amin but 
still a potential barrieio near a = just as in Sect. (2), and choose t] such 
that Veff is real and positive near amax- An example is obtained by choosing 
T] = 0.7, ( = 0.02 which gives an amax near a = 8.1. Fig. (7) shows Veff near 
amax, showing the positive barrier as one approaches from the left. 



'^We are assuming that the approximation (l26l) holds also near a = 
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In summary, the points amin and ttmax in the quantum equation imply 
turning points in the effective classical dynamics (126 p in those cases where 
Veff remain real. Thus the modified Friedmann equation can support a 
cyclic universe; for a review on cyclic universes see [13]. Since the difference 
equation f lT9|) does not depend on the sign of a and hence represents both 
expanding and contracting solutions, one may argue that the underlying 
quantum universe is also cyclic. 

Finally, we may also use the effective classical dynamics to interpret the 
oscillations seen in probability density curves: they imply that the expansion 
of the classical FRW-A universe will deviate from the pure exponential, a{t) ~ 
exp(t(l — S)), with 6 the small oscillatory term discussed in the last section. 

7 A = Case 

In a fiat universe without a cosmological constant (Minkowski spacetime) the 
linear WDW is just the free time- independent Schrodinger equation whose 
solution is ip{a) = cia + c^ where c, are complex constants. The probability 
density is then j9(a) = Aa^ — Ba + C with A, C > 0. For i? < the solution 
is not normalisable but for sufficiently large i? > one gets p < beyond a 
certain a = amax and hence an allowed region < a < amax- 

When A = 0, the a^ piece is absent in the nonlinear difference equation 
fin?]) and we set the constant a = — 1 in flTTl) . As preliminary numerical 
trials indicated that p increased without bound in the forward evolution if 
the initial p were nonvanishing, we report the results when we set p^ ~ at 
a = 10 or a = 20 and choose p to be a small positive value ~ 0.2. In Fig. (8) 
we see also that the nonlinearity makes the p curves less steep compared to 
that of the linear WDW equation for the same initial conditions at a^nax ■ 

For 7] = 0.2, 0.5 there is no potential barrier formed for small C, but just 
as in the case for A 7^ there exists a Cc beyond which the quantum universe 
has a minimum allowable size, in addition to a maximum size at a = 10 or 
a = 20 that we had set by hand, see Fig. (9). 

In passing we note that the exact solution fl2^ does not exist for the 
case A = because when the a^ piece is absent the equation is inconsistent. 
However this does not mean that solutions with both an amin > and a^nax 
do not exist when A = but rather that one must either consider more than 
three lattice points, or realise that for a given a^ and z the actual solution 
might have p vanishing before the chosen p_|_ = or after the chosen p_ = 0, 
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as the numerical results do indicate, Fig. (9). In other words a three lattice 
point solution may be constructed, for example, by first choosing p_ < and 
then restricting the solution to the region where the probability density is 
positive. 

8 Conclusion 

The probability density describing a quantum de Sitter universe in the non- 
linear WDW framework of Ref.jl] was found to obey a nonlinear difference 
equation. For an arbitrarily weak nonlinearity, (, the universe has a maxi- 
mum allowable size, amax ~ ^, while a potential barrier screens the classical 

a = region when the free parameter rj < 3/4. 

As ( increases beyond Cc? which is still very small, there emerges a minimal 
allowable size Omm to the quantum universe for any value of t]. One may 
think of the quantum nonlinearity as counter-acting [9] the dispersion of the 
wavefunction, localising it to a finite range amin < a < cimax- 

In the quantum picture the existence of amin and a^ax was implied by 
the vanishing of the probability density: The consistency of interpreting amin 
and ttmax as such was seen by examining the effective classical dynamics. We 
showed that near such nodes of the probability density there correspond 
bounces in the effective classical evolution for cases where Vg// remains real; 
in some cases we have a self-consistent interpretation of a cyclic universe. 

While a modified dynamics replacing the Big Bang by a Big Bounce might 
have been expected, it is remarkable that the evolution of the universe at large 
scales has also been modified by the nominally weak quantum nonlinearity 
F: Instead of an unending exponential expansion, there can now be a Big 
Crunch. This significant modification of the classical dynamics is due to 
the enhancement ^ of the nonlinearity F near the nodes of the probability 
density. 

The current accelerating phase of our universe is well modeled by a domi- 
nant cosmo logical constant. We interpret our results as suggesting that such 
an accelerating phase cannot continue indefinitely but will be replaced by a 
contracting phase. Since our information-theoretic approach does not pre- 
sume a specific micro-dynamics responsible for the nonlinearity F, it might 
be summarising various possibilities. For example, in Ref.[T5] a semi-classical 
analysis showed that in gravitational theories with extra dimensions, such as 
those that appear in M-theory frU\, a vacuum state with non-zero cosmologi- 
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cal constant is generically unstable, one possibility being an eventual collapse 
in a Big Crunch: Our results are compatible with those of ^5\. In this re- 
gard, it is interesting to note that cyclic models of the universe [17] inspired 
by M-theory y^ have been proposed. 

Of course to truly see a cyclic time-evolution in the quantum domain of 
our model we need to go beyond the FRW-A case and include, for example, a 
massless scalar field to act as an internal clock [12], the FRW-0 model. The 
study of the FRW-0 model is important since the classical Big Bang has a 
curvature singularity. 

Based on our perturbative study of the latter model in Ref. [1], and our 
analysis of the nonlinearity F in Section(5), it is plausible that a non-singular 
cyclic evolution of the quantum universe might also be seen for the FRW-^ 
model within our nonlinear WDW framework [TT]. Also, a cyclic evolution 
should also hold for models in which the cosmological "constant" is slowly 
varying, in particular for some realistic inflationary potentials [TT] . 

Ultimately the free parameters rj and ( need to be fixed by comparing 
computations of observables with empirical data. On the other side, the 
robustness of the theoretical conclusions may be checked by examining in- 
formation measures which are deformations of the KL measure ^7j| . 

Finally we note that the nonlinear difference equation ( !T9|l may be viewed 
as a one-dimensional nonlinear map and could possibly exhibit chaos for 
certain range of parameters [20] • It might be useful to investigate the physical 
relevance of such behaviour for the cosmology considered. 
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Figure Captions 

• Figure 1 : The barrier in the effective potential for ( = 5 x 10~^ < (c- 
(The region below a = 5 x 10"'^ has been extrapolated since Q cannot 
be evaluated there using the central difference formula). 

• Figure 2: Probability density curve for ^ = 5 x 10^'^. For low a the 
curve follows very closely the curve for the linear theory but p becomes 
complex at a = 14.4. 

• Figure 3 : The region of Fig. (2) near ttmax- 

• Figure 4 : Variation of a^ax with ( for t] = 1/2. 

• Figure 5 : The probability density curve for ( = 0.1 > (^ = 0.049 
showing both an Omm = 4.55 and amax = 5.3. The discrete data points 
have been connected for visualisation. Included are the two unphysical 
points at the ends. 

• Figure 6 : The oscillations in the p curve shown here in the ratio plot 
of p/pasymp whcre Pasymp IS the asymptotic value for the linear theory, 
Eq.®. 

• Figure 7 : Veff near ttmax for rj = 0.7, ( = 0.02. The last point on the 
right has imaginary components and is unphysical. 

• Figure 8: Curves for A = 0. The upper curve is for the linear WDW 
equation while the lower is for the nonlinear WDW equation with the 
same initial conditions near a = 20. 

• Figure 9: An example for A = showing both an amin and ttmax- Note 
the unphysical points at the ends. 
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